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Abstract 

The aim of this study was to develop a novel technique for acquiring 3-dimensional (3D) coronary CINE magnetic resonance 
images with both water and fat visualization during free breathing and without external respiratory or cardiac gating. The 
implemented multi-echo hybrid 3D radial balanced Steady-State Free Precession (SSFP) sequence has an efficient data 
acquisition and is robust against motion. The k-space center along the slice encoding direction was repeatedly acquired to 
derive both respiratory and cardiac self-gating signals without an increase in scan time, enabling a free-breathing 
acquisition. The multi-echo acquisition allowed image reconstruction with water-fat separation, providing improved 
visualization of the coronary artery lumen. Ten healthy subjects were imaged successfully at 1.5 T, achieving a spatial 
resolution of 1.0x1.0x3.0 mm^ and scan time of about 5 minutes. The proposed imaging technique provided coronary 
vessel depiction comparable to that obtained with conventional breath-hold imaging and navigator gated free-breathing 
imaging. 
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introduction 

Coronary magnetic resonance angiography (CMRA) iias been 
highly desired because of its potential clinical significance for 
managing coronary artery disease [1]. However, the features of the 
coronary arteries, such as the tortuous shape, the small vessel size, 
and the sensitivity to motion make routine CMRA still challenging 
in current clinical practice. To prevent motion artifacts, conven- 
tional CMRA acquisitions are usually triggered by the electrocar- 
diographic (ECO) signal and data acquisition is confined to the 
mid-diastolic period of the cardiac cycle when the heart is 
relatively at rest. As this period is variable with a changing heart 
rate, a precise subject-dependent ECG trigger delay needs to be 
determined. Different segments of the coronary anatomy may also 
be quiescent at different points in the cardiac cycle [2]. To 
overcome these problems, the ability to image the coronary 
arteries throughout the cardiac cycle (coronary CINE imaging) is 
highly desired. 

Effective suppression of the epicardial fat signal is important in 
CMRA to improve coronary vessel conspicuity [3]. In conven- 
tional CMRA, a long spectrally selective RF pulse centered on the 
fat resonance frequency followed by spoiler gradients is typically 
used for this purpose. The susceptibility of lung-heart interface can 
significantly decrease the field homogeneity, causing incomplete 
fat suppression. Furthermore, such approach interrupts the 
continuous data acquisition during steady-state, making CINE 
acquisition difficult. Several alternative fat suppression techniques 
have been applied to achieve coronary CINE imaging, such as 



contrast-enhanced or phase-detection techniques [4—7]. Effective 
fat suppression pulses have been apphed to coronary CINE 
imaging by mitigating the steady-state interrupts [8]. Linear 
combination SSFP (LC-SSFP) method was shown to achieve 
simultaneous imaging of the coronary lumen and the epicardial fat 
which is advantageous in enhancing vessel conspicuity at the cost 
of longer acquisition time [9] . 

In this study, a novel free-breathing multi-contrast coronary 
CINE imaging technique is described that provides coronary 
water and fat CINE images efificiendy without requiring dedicated 
fat suppression pulses. This technique is based on a multiple-echo 
hybrid 3D radial balanced steady-state free-precession (SSFP) 
sequence [10], which has been demonstrated to be a promising 
technique for achieving self-gated free-breathing 3D cardiac CINE 
imaging for functional measurement and further exploited with 
water and fat separation in this study. The performance of the 
proposed technique is evaluated by comparing with the conven- 
tional breath-hold CMRA and free-breathing navigator gated 
imaging techniques in healthy volunteers. 

Methods 

This study was approved by the WeiU-Cornell Medical Center 
institutional review board and performed in compliance with the 
Health Insurance Portability and Accountability Act (HIPAA). 
Written consent was obtained from all subjects prior to imaging. 
AU imaging experiments were performed on a 1.5T GE HDx 
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scanner (General Electric Healthcare Technologies, Waukesha, 
WI, USA). 

Self-gated Multi-echo 3D CINE Radial SSFP Sequence 

We implemented a multi-echo stack-of-stars 3D balanced SSFP 
coronary CINE sequence with conventional Cartesian sampling in 
the kj (slice encoding) direction and radial sampling in the k^^-ky 
plane similar to that developed previously for cardiac functional 
imaging [10]. A unique feature of this sequence is that the 
acquired k-space center line along the k^ direction can provide a 
direct measure of the global motion of the imaged volume cdong 
this direction, from which the superimposed respiratory and 
cardiac motion can be separated by filtering and used to perform 
retrospective motion gating of the acquired image data [10]. 
Unlike conventional ECG-gated coronary imaging sequences, 
which only capture the coronary arteries in a single cardiac phase 
(typically during mid-diastole), the proposed sequence acquires 
data continuously throughout the cardiac cycle and is capable of 
providing multiple cardiac phase resolved images, from which an 
image with the optimal depiction of the targeted coronary artery 
could be selected. A major technical challenge in such continuous 
data acquisition scheme is effective epicardial fat suppression 
which is vital for optimal coronary vessel depiction. In this study, 
the multi-echo data acquisition was exploited to provide simulta- 
neous water and fat images of the coronary arteries, obviating the 
need for a separate fat suppression pulse. 

Water-fat Separation for Multi-echo 3D CINE Radial CMRA 

The multi-echo radial acquisition is capable of separating the 
water and fat signals by linearly combining data sets acquired at 
different echo times [11-16]. Specifically, the MR signal </, 
obtained from a single voxel in the «-th echo data set was modeled 
as follows (ignoring the local static field inhomogeneity) [12]: 

di = w + (j^''^-='^'^i-f (1) 

where w and f denote the water and fat signals respectively, TE, is 
the f-th echo time, and fc, denotes the known chemical shift 
between water and fat peaks (about —210 Hz at 1.5T). The 
unknowns w and / can be obtained by inverting the following set of 
Unear equations with a 3-echo acquisition: 



1 ^^^fcsTE^\ /dl\ 
1 e/2-fe7-£2 L2 (2) 

In our study, TEs were chosen such that TE2~TE] = TEg~TE2 
(TEs were 0.3, 1.6, 2.9 ms, respectively). 

Non-linear Combination of Water and Fat Images to 
Selectively Improve Coronary Lumen Visualization 

As described above, we can simultaneously obtain water and fat 
images of the coronary artery. The water-fat separation usually 
suppresses the main spectral peak of the fat signal (e.g., —210 Hz 
at 1.5T) in the water image. However, the other smaller fat 
spectral peaks may not be well suppressed, which could degrade 
the contrast between the coronary artery and the surrounding fat 
[17]. In this study we proposed a simple and fast non-linear 
combination method of water and fat images to achieve improved 
coronary lumen imaging as follows: 



W,u„en = H-(\ (3) 

V max(\f\)) 

where w and / are the separated water and fat images respectively, 
|-| denotes the absolute value (image magnitude), and max(\f\) is 
the manually selected regional maximum intensity of the 
epicardial fat surrounding the coronary artery. The non-linear 
image intensity transformation in Eq. 3 scales the water image 
pixel-wise with a weighting mask derived from the fat image, 
which has the effect of further suppressing the residual fat signals 
in mostly fat voxels in the water image where |/| wwax(|/|). 

Phantom Scans 

To demonstrate the water-fat separation and combination 
techniques in vitro, a water-fat phantom consisting of a balloon 
filled with peanut oil and immersed into a water container was 
imaged using the developed multi-echo 3D radial balanced SSFP 
sequence with the following imaging parameters: TR = 4.9 ms, 
three TEs = 0.3/ 1 .6/2.9 ms, FA = 40°, BW = ± 1 25 kHz, 
FOV = 29cm, spatial resolution = 1.1 x 1.1 x5.0 mm^, 16 shoes, 
400 in-plane projections, scan time of 31 s. The volume body coil 
was used for signal reception. 

Human Scans 

3D CMRA was performed in 10 healthy volunteers (mean 

age = 29.4 years ±4.2, 9 male, 1 female). An 8-channel cardiac 
phased array was used for signal reception. The right coronary 
artery (RCA) and the left anterior descending coronar)' artery 
(LAD) were imaged in targeted double-oblique planes in 
randomized order. The imaging parameters of the self-gated 
(SO) multi-echo 3D radial balanced SSFP CINE CMRA sequence 
were as follows: TR = 5ms, three TEs = 0.3/1.6/2.9 ms, 
FA = 40°, BW = ± 1 25 kHz, FOV = 26 cm, spatial resolu- 
tion= 1.0x1.0x3.0 mm^, 14-16 slices, free-breathing acquisition 
without ECG gating. The- temporal resolution (the length of each 
cardiac phase) was 70-80 ms by sharing small amount of high 
frequency data from adjacent cardiac phases [10,18]. The scan 
time was approximately 5 minutes. Retrospective respiratory 
gating was used to reconstruct CMRA images with a 50% self- 
gating efficiency. 

For comparison, 3D CMRA was also obtained with the 
conventional breath-hold (BH) ECG-triggered 3D Cartesian 
balanced SSFP sequence as well as a free-breathing diaphragm 
navigator (NAV) gated ECG-triggered 3D Cartesian balanced 
SSFP sequence [19]. The BH sequence was applied with the 
following imaging parameters: TR = 4.2 ms, TE= 1.7 ms (partial 
echo), FA = 60°, BW= ±62.5 kHz, FOV = 26 cm, spatial resolu- 
tion = 1.0x1.6x3.0 mm^, 8 slices, 200 ms data acquisition win- 
dow during the mid-diastolic coronary rest period. The scan time 
was 24 heartbeats. A subject-specific ECG trigger delay from the 
R wave to the start of the data acquisition was determined from 
2D CINE scout images. The NAV sequence was apphed witii the 
following scanning parameters: TR = 4. 1 ms, TE = 1.6 ms (partial 
echo), FA = 60°, BW=±64kHz, FOV =26 cm, spatial resolu- 
tion= 1.0x1.0x3.0 mm'^, 14 slices, 128 ms mid-diastolic data 
acquisition window, same ECG trigger delay as the BH CMRA 
scan. The phase ordering with automatic window selection [20] 
real-time navigator gating algorithm was used with a 5 mm gating 
window. Spectrally selective RF pulses were used for epicardial fat 
suppression for both BH and NAV CMRA sequences. The three 
CMRA sequences (SG, BH, and NAV) were performed in a 
random order for all the subjects. 
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Data Analysis 

For the water-fat phantom, region-of-interest (ROI) analysis was 
used to measure the contrast between the water and fat signals. 
For in vivo CMRA, the contrast between coronary artery (RCA/ 
LAD) and surrounding fat was measured. The coronary ROI was 
placed on the proximal segment of each coronary artery branch 
and the fat ROI was placed on the surrounding epicardial fat. The 



contrast was calculated as c = 



where Iw and L- are the 



average water and fat signal intensities within the ROIs. The 
change in contrast by applying the water and fat image 
combination (Eq. 3) was calculated as ^^f^ x 100%, where Ci 
and 0-2 are the measurements without and with the combination, 
respectively. In this study, we selected the best phase of the cine 
SG data based on the overall subjective vessel visualization (which 
considers factors such as coronary sharpness, length, residual 
motion artifacts and contrast). This optimal cine SG image was 
used for the measurements in all tables. 

Images were reformatted using OsiriX imaging software 
(OsiriX Foundation, Geneva, Switzerland). For each imaged 
coronary artery, an experienced radiologist scored the vessel 
depiction using a 5-point scale: 4-exceUent, 3-good, 2-fair, 1-poor, 
0-not seen. Wilcoxon paired sample signed-rank test was 
performed on the scores to assess the difference between the 
CMRA sequences. A significance level of less than 0.05 was 
considered to indicate statistical significance. To evaluate the 
vessel sharpness, the local maximum (/^J and minimum (/„„„) 
intensity values across the each side of the vessel were determined, 
from which image sharpness was calculated as the inverse of the 
distance between 0.8f4„;,-4i„j+/„i„ and 0.2C4„-4,„J+4m [21,22]. 
A cross-section of the vessel at the proximal RCA/LAD was 
chosen for measuring the vessel sharpness. The sharpness obtained 
at two sides was then averaged for evaluation and comparison. 
When the vessel sharpness cannot be measured due to excessive 
motion artifacts, a sharpness of 0 will be assigned. The vessel 
length of depiction was also measured and evaluated. 

The start time, duration and best phase of the rest period were 
measured as percentages of the cardiac cycle duration for different 
subjects and coronary branches. Statistical analysis on the 
measurements (mean and standard derivation) was reported. 
The measurement difference (p-value) between RCA and LAD 
was also reported. 

Results 

Figure lA-C shows the separated water and fat images as well 
as the combined image with Eq. 3 of the water-fat phantom. The 
combined image (Figure IC) demonstrates a better contrast (48% 
improvement) than those in the separated water image (Figure 1 A). 
Figure ID displays a signal cross-section of the water and 
combined images (Figure lA&C), demonstrating the image 
intensity changes by applying the water-fat combination (Eq. 3). 

CMRA was obtained successfully in all 10 subjects. Figure 2 
shows an example of water and fat CINE images of the RCA using 
the proposed technique at five representative cardiac phases 
centered around the mid-diastole, demonstrating excellent coro- 
nary artery delineation and effective cardiac and respiratory 
motion suppression in both water and fat images. 

In Figure 3, water, fat and the combined image with enhanced 
coronary lumen (Eq. 3) are shown for both RCA and LAD images. 
Note the improved contrast of the coronary artery lumen against 
the surrounding fat (Figure 3C). 

Table 1 summarizes the coronary artery contrast of RCA and 
LAD averaged over 10 subjects, demonstrating significandy 
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Figure 1. Phantom water image, fat image, and combined 
image obtained with Eq. 3 are shown in A-C. The central 
horizontal lines of A and C are plotted in D. 
dol:10.1371/journal.pone.0089315.g001 



improved coronary contrast after non-linear combination of water 
and fat images (Eq. 3), compared to water-only images. 

Comparisons of BH, free-breathing NAV, and our proposed SG 
free-breathing techniques are demonstrated in Figure 4. The SG 
method with water-fat separation generated both water and fat 
images (only water images are showed as the third column in 
Figure 4). The NAV scan time was 5.1±1.3min and the 
diaphragm navigator gating efficiency was 43 ± 1 2 % . The heart 
rate was 60.0±6.5 bpm. Both right and left coronary arteries are 
well depicted. SG method was found to provide slighdy lower but 
statistically comparable (not statistically different) vessel depiction 
scores when compared to the conventional BH and free-breathing 
NAV techniques (Table 2). 

For the coronary CINE images obtained with SG method, we 
have measured vessel sharpness of both RCA and LAD during all 
cardiac phases where the vessels were not corrupted by motion 
(referred as rest period). The sharpness values at cardiac phases 
where the vessel was corrupted by motion were assigned with 
zeros. Figure 5 shows the change in coronary vessel sharpness 
throughout the cardiac cycle for all 10 subjects. Note that the 
cardiac phases start at end-systole. The start time, duration and 
best phase of the rest period vary for different subjects and 
coronaries. The statistical analysis of the vessel rest period is 
summarized in Table 3. It shows that the LAD rest period starts 
earlier and lasts longer than that of RCA. RCA rest period has a 
large variation among subjects. 

Vessel sharpness and length of depiction were measured for BH, 
NAV and SG methods. For SG method, only the best cardiac 
phase was analyzed. The results are summarized in Table 4, 
showing similar performance among the three techniques (all P 
values > 0.05, N= 10). 

Discussion and Conclusion 

It is well known that there exists a large variability in starting 
point and duration of the coronary rest period (depending on 
subject and coronary branches) with no correlation to heart rate 
[2]. Conventional single-cardiac-phase coronary imaging would 
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Figure 2. Respiratory and cardiac self-gated free-breathing coronary water (top row) and fat (bottom row) cine images are 
displayed at five representative phases (every other phase). The right coronary artery Is clearly seen In mid-dlastole (dashed box) In both 
water and fat Images, and visualization Is sensitive to cardiac phase. 
dol:10.1371/journal.pone.0089315.g002 



ideally require the precise estimation of patient- and vessel branch- 
specific trigger delay and acquisition window, which is difficult and 
can be time-consuming. Our preliminary results have demon- 



B 



strated the feasibility of simultaneous water and fat 3D CINE 
imaging of the coronary arteries during free breathing without the 
need for external motion gating devices. Given the segment-to- 
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Figure 3. Multi-contrast coronary imaging (top row: RCA, bottom row: LAD). Water Images are shown In A&D, fat Images In B&E, combined 
Images with enhanced coronary artery contrast (Eq. 3) In C&F. 
dol:1 0.1 371 /journal.pone.008931 5.g003 
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Table 1. Coronary artery to epicardial fat contrast 
measurements and their corresponding changes without and 
with the non-linear water and fat image combination (Eq. 3) 
(N = 10). 







Artery to fat 








contrast 


Water image 


Combined image 


Change 


RCA 


0.7 ±0.1 


1.1 ±0.3 


*52.1% ±38.4% 


LAD 


0.6 ±0.7 


1.0±0.6 


•42.6% ±53.0% 



Statistically significant differences (P<0.05} are marked with*. 
doi:l 0.1 371/journal.pone.008931 S.tOOl 



segment variation in coronary motion through tlie cardiac cycle, 
CINE imaging permits optimal retrospective selection of the 
cardiac phase that provides the best visualization for any given 
coronary segment, a major benefit over existing single cardiac 
phase CMRA techniques. Future studies in a larger population 
wiU be needed to demonstrate this potential advantage of the cine 
imaging approach. 

The self-gating method apphed in this study has been vahdated 
for cardiac CINE imaging for functional measurement, by 
comparing to the standard external measures [10,23]. Although 
the orientation of the targeted coronary artery imaging might be 
dilferent from that of the cardiac CINE imaging, we have 
observed reliable gating signals in our coronary imaging studies. 
Respiratory motion of the heart is mainly along the head-foot 
direction. Although the motion projecting to direction of an 
arbitrary oblique plane may have an under-estimated amplitude, 
the respiratory motion waveforms (motion frequencies) remain. 
Thus the self-gating signal based on the center-of-mass can still be 
derived with a band-pass filtering. In this study, a fixed respiratory 



Table 2. Comparison of coronary vessel depiction obtained 
with breath-hold (BH), free-breathing NAV (NAV), and self- 
gated (SG) free-breathing CMRA techniques (N = 10, P>0.05). 







Vessel depiction 


BH 


NAV 


SG 


RCA 


3.6±0.7 


3.1 ±1.4 


3.0±0.7 


LAD 


3.4±1.0 


3.2±1.1 


2.8±0.9 





doi:10.1371/journal.pone.0089315.t002 



gating efficiency 50% was applied [10], however a subject-specific 
gating efficiency needs to be further investigated. 

In this study, the coronary image quality was found to be similar 
for the proposed SG technique and the conventional BH and 
NAV techniques, although these sequences had different temporal 
resolutions of 78 ms, 200 ms and 128 ms, respectively. This may 
be explained by the relatively lower heart rate (correspondingly 
longer mid-diastolic cardiac rest period) and the ability to sustain a 
long breath-hold (24 heartbeats) in young healthy volunteers of 
this study. Comparison of image quality in a patient population 
will be considered in our future work. 

Arrhythmia was not obseived in the healthy volunteers 
recruited for this study. Heart rate has a direct effect on the 
acquisition time of the breath-hold sequence and together with 
respiratory motion influences the acquisition time of the 
diaphragm navigator gated sequence. The acquisition time of 
the BH sequence changed according to the heart rate, since a fixed 
number (24) of heart beats were acquired. We did not observe a 
significant correlation between the heart rate and scan time for the 
navigator-gated sequence. The acquisition time of the self-gated 
radial sequence was fixed at 5 minutes. 

The quality of the coronary MRA is primarily determined by 
the effectiveness of motion suppression (both cardiac and 
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Figure 5. Vessel sharpness obtained with SG technique was measured throughout the cardiac phases. A sharpness of 0 indicates that 
the vessel sharpness cannot be measured due to excessive motion artifacts. 
dol:1 0.1 371 /journal.pone.008931 5.g005 



respiratory) and to a lesser degree by the spatial resolution. The 
higher quality of the breath-hold sequence as observed in this 
study, while surprising in light of previous studies [24], can be 
explained by our special study population comprising of young, 
healthy and seasoned volunteers who have low heart rates (hence 
minimizing cardiac motion effects) and are capable of sustaining a 
long breath-holds of 24 heartbeats. The images were reviewed by a 
radiologist with five years experience in cardiovascular MRI. As 
our study population consists of young seasoned volunteers 
capable of sustaining good-quality long breath-holds, the reader 
gave higher scores to breath-hold images which generally have less 
respiratory motion artifacts when compared to the navigator gated 
images. We expect a different result in an older patient population, 
which win be studied in subsequent works. 

The current temporal resolution (~ 70-80 ms) with the 
proposed SG technique is not high enough to provide nice 
coronary images throughout all the cardiac phases (Figures 2&5). 



Table 3. Statistical analysis of RCA and LAD rest period. 





% of the 
cardiac cycle 


Best Phase 
location 


Start Time 


Rest Duration 


RCA 


54.7 ±19.0 


44.1 ±6.4 


35.9±7.2 


LAD 


48.6±9.5 


40.3±6.2 


39.6±7.2 


RCA v.s. LAD 


P = 0.19 


*P = 0.02 


P = 0.38 


All numbers are report as percentages of the cardiac eye 
difference between the measurements of RCA and LAD is 


e duration. The 
marl<ed with *, if the 



difference is statistically significant {P<0.05). 
doi:10.1371/journal.pone.0089315.t003 



But with radial acquisition and view sharing reconstruction, the 
SG technique could provide reasonable coronary images at 
multiple cardiac phases rather than one or two using conventional 
methods (BH or NAV). Since the self-gating method is based on 
the blood intensity change during the cardiac cycle, the cardiac 
phases showed in Figure 5 start at end-systole (with the smallest 
amount of blood signal) [10]. Table 3 shows that the LAD rest 
period starts earlier and lasts longer than that of RCA, which 
agrees with what was reported earlier [2]. The large variance of 
RCA rest period (±19% of the cardiac cycle) among different 
subjects indicates that patient-dependent trigger delay is crucial for 
the conventional approach using an ECG-triggered single-phase 
imaging sequence, while it is much less problematic for the multi- 
phase cine sequence developed in this work. 

Table 4. Comparison of coronary vessel length of depiction 
and sharpness obtained with breath-hold (BH), free-breathing 
NAV (NAV), and self-gated (SG) free-breathing CMRA 
techniques {N = 10). 







Measurement 




BH 


NAV 


SG 


Length (cm) 


RCA 


12.6±2.2 


13.2±1.6 


11.5±2.3 




LAD 


10.5±2.5 


9.6±2.5 


10.2±1.3 


Sharpness 


RCA 


0.92±0.26 


1.00±0.10 


0.93 ±0.1 5 


(mm"') 












LAD 


0.87±0.15 


1.05±0.24 


1.01±0.14 





The p-values between methods are great than 0.05 (the difference is not 

statistically significant). 

doi:1 0.1 371 /journal.pone.008931 5.t004 
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A simple model of water-fat separation without considering the 
local field inhomogeneities was used in this study. When used at 
1.5 T field strength in combination with targeted volume 
shimming to reduce field inhomogeneities in the heart region, 
our technique was shown to consistently provide robust epicardial 
fat suppression and coronar}' angiograms of good quality in our 
study cohort comprised of young healthy volunteers. However, 
water-fat separation with simple Dixon method can be problem- 
atic when good shimming is difficult to achieve (e.g., in larger 
patients), and a more advanced Dixon method may be needed. 
The comparison of our technique and more advanced water-fat 
separation methods that take into account local inhomogeneity, 
such as in generalized three-point Dixon, IDEAL, or VARPRO 
methods [12,13,16], will be considered in our future work. 

As reported in Tables 1&4, the decreased values of the SG 
method are not statistically significant compared to the other two 
methods. However, we have to note that selecting the best rest 
period of the cardiac cycle with the proposed method may not lead 
to the best image quality, which may also be affected by many 
other aspects, such as respiratory motion, sampling trajectory, data 
undersampling, fat-sat strategy, and so on. Although a fairly large 
50% gating window was intentionally chosen in this study to keep 
the acquisition time for SG and NAV methods similar, a narrower 
gating window (20~30%) may be used to improve coronary 
quality if more scan time is available. The Dixon method used in 
our study was based on multi-echo radial acquisition. Although we 
had corrections for gradient delays and eddy currents, the residual 
errors may still affect the radial trajectories within a single echo 
data set or among different echo data sets and result in image 
quality degradation. Thus, a different fat-sat strategy used with 
respect to BH and NAV could also influence the vessel sharpness 
and the subjective score. 

In this study, a nonlinear filtering was used to enhance the 
contrast of the coronary artery by combining the simultaneously 
obtained water and fat images (Eq. 3). A region of epicardial fat 
surrounding the coronary artery' was manually selected and its 
averaged intensity was calculated as the regional maximum 
intensity for Eq. 3. Although this value is just an estimation of 
the fat signal intensities around the vessel, this nonlinear filtering 
method is relatively insensitive to the inhomogeneity of the fat 
signal, especially compared to a simple linear mask based on a 
threshold of the fat signal (data not shown). Because of the 
nonlinear filtering, the background noise of the combined images 
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